Destruction of the small Fermi surfaces in Na3;Co02 by Na disorder 
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We show using density functional calculations that the small e'g Fermi surfaces in Naa:Co02 are 
destroyed by Na disorder. This provides a means to resolve the prediction of these sections in band 
structure calculations with their non-observation in angle resolved photoemission experiments. 
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The layered oxide, Nai^Co02, has attracted consid- 
erable interest because of its enhanced thermoelectric 
properties^ and because, when hydrated the material 
becomes a superconductor, possibly with unconventional 
pairing:^'-. Naj:Co02 consists of triangular sheets of Co 
ions, with nominal d electron count, 5-t-a;. These are 
coordinated by distorted edge sharing octahedra, made 
from triangular sheets of O ions above and below the Co 
sheets. These C0O2 tri-layers are stacked in an alter- 
nating fashion to form prismatic sites, which contain the 
Na. 

Establishing the electronic structure is requisite for un- 
derstanding of both the thermoelectric behavior and the 
superconductivity. Local density approximation (LDA) 
band structure calculations, ''■^ show a Fermi level that 
lies near the top of a narrow manifold of t2g bands, which 
is separated by a gap from the higher lying tg derived 
states. With the actual rhombohedral Co site symmetry, 
the three t2g orbitals are further divided into an ag and 
two degenerate Cg symmetry orbitals, denoted e' to dis- 
tinguish them from the higher lying Cg manifold defined 
by the primary crystal field splitting. Significantly, local 
density approximation band structure calculations show 
the presence of two sheets of Fermi surface over a wide 
range of x. These are a large ag derived hole cylinder 
around T and six small hole sections along the hexagonal 
T - K lines. The calculated Fermi surfaces of the hy- 
drated superconducting compound similarly show both 
ttg and small e' sectionsi^ While these e' sections are 
small, and therefore are not expected to contribute sub- 
stantially to the conduction, they contribute strongly to 
the density of states and should be very important for 
scattering, magnetic susceptibility and other properties 
because of the two dimensional nature of the material 
and their heavy mass.^'* 

In contrast to the band structure calculations, angle 
resolved photoemission (ARPES) experiments by sev- 
eral groups at various doping levels observe the large Og 
cylinder but do not find the small e' sectionsAiSiiLiSii^ 
Early on,^ it was noted that although metallic, Na2:Co02 
had a small band width compared to plausible values 
for the on-site Coulomb repulsion, U . One explanation 
that has been advanced for the absence of the e' sec- 
tions is that they are destroyed by correlation effects. 
LDA-I-U calculations, which incorporate Coulomb repul- 
sion in a static mean field like way, have indeed shown 



that the small sections are removed for reasonable values 
of [/ii4 However, the LDA-I-U approximation also fa- 
vors insulating, charge ordered ground statesi Removal 
of the small sections is also found in the large U limit 
Gutzwiller approximation^^ However, in metals, fluctu- 
ations, including charge and orbital fluctuations are fa- 
vored by kinetic energy considerations and the presence 
of e' sections would open more degrees of freedom for 
fluctuations. Dynamical mean flcld calculations, which 
incorporate physics like the LDA-I-U approximation, but 
in a beyond mean fleld manner that incorporates fluc- 
tuations, yielded the opposite tendency to the LDA-I-U 
approximation, i.e an enlargement of the e' pockets over 
the LDA."'^^ Low temperature speciflc heat measurements 
on superconducting hydrated Naa;Co02 show two dis- 
tinct energy gaps,"'^ which is difficult to understand if 
there is only one, simple Fermi surface. Thus, the dis- 
crepancy between the ARPES and LDA Fermi surfaces 
remains a puzzle. Here we resolve this puzzle by showing 
that the small sections are destroyed by scattering due 
to the Na disorder inherent in the structural chemistry 
of Na:rCo02. 

As mentioned, the Na ions are contained in O prisms 
defined by the O triangular lattice. There are two types 
of sites: Those directly above a Co ion, denoted Naco 
here, and those above the holes in the Co sheet, denoted 
Na^, (these are directly above the O atoms in the O layer 
on the opposite side of the Co sheet). There is one Naco 
and one Na/i site per Co, with a total Na filling of x per 
Co, or x/2 per site. However, the number of available 
sites is smaller since occupation of an Naco next to an 
occupied Na^ site or vice versa would be highly unfavor- 
able due to the short Na - Na distance of 1.64 A that 
would result. In fact, local structures that within a given 
Na layer have Na of a given type (Naco or Na/j) hav- 
ing Na neighbors of the same type are favored. Also Na 
in adjacent layers tend not to coordinate the same Co 
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To investigate the effect of Na ion ordering, we focus on 
the X = 2/3 composition and perform electronic structure 
calculations for two supercells with different Na order- 
ings and Co sites. Specifically, we report results for two 
\/3 X \/3 cells with lattice parameters a=2.84 A, c=10.81 
A. Results are reported for the ideal structure, with api- 
cal O height zo=0.0864 (Ref. H), but with different Na 
orderings, and for fully relaxed atomic positions. Each 




FIG. 1; (color online) Structure of two supercells, A (left) and 
B (right). Na ions are indicated by the large blue spheres, Co 
by the mid-size green spheres, and O by the small red spheres. 
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FIG. 2: (color online) Band dispersions along the T-Ad and 
A-L directions for the two relaxed supercells: A (top) and B 
(bottom). The different line types are for the three V-M or 
A-L directions, separated by 120°, which are non-equivalent 
in the relaxed supercells. 



supercell contains 6 Co and 12 Na sites in two layers. 
The two supercells, as shown in Fig. nare: (A) with one 
Na layer on the Naco sites and the other on the Na/j sites 
(both with two Na atoms per layer); and (B) with both 
Na layers on the Na;i sites. 

The calculations were done within the local den- 
sity approximation with the general potential linearized 
augmented planewave (LAPW) method including local 
orbitalsiSiSSiS^ Well converged basis sets, of approxi- 
mately 2600 LAPW functions plus local orbitals, were 
used, with LAPW sphere radii of 1.96 aq for the Na and 
Co and 1.52 ag for O. 

The calculated band structures near the Fermi energy, 
Ep, for the relaxed cells are shown in Fig. [3 along the 
supercell T-M lines. Due to the -s/S x \/3 supercell pe- 
riodicity, these are folded T-K directions in the original 
hexagonal zone, and therefore intersect the small pock- 



ets. The band structures obtained for the two supercells 
are very similar to each other and to the virtual crystal 
band structurei^ In particular, at this doping level there 
is sufficient k^ dispersion that only the outer of the two 
flg and e' pairs of bands produce Fermi surfaces in the 
kz = plane, while both ag and neither e' bands pro- 
duce Fermi surfaces at kz = 1/2. This differs in detail 
from virtual crystal calculations^^ where both e' bands 
reach maxima above Ep in the k^ = 1/2 plane. This 
difference in the e' , but not the Ug topology, when using 
real Na ions instead of the virtual crystal approximation 
is an indication that variations in the Na potential due 
to disorder could localize the e'g carriers. In any case, 
the maximum of the Ug band is 0.17 eV above Ep, while 
the e' maximum is only 0.01 eV above Ep. Besides the 
positions of the band maxima, there is also a large differ- 
ence in band velocity along the T-M line, i;(ag)=1.4 eV 

A, v{e'g)=0.3 eV A. 

The second indication of localization of the e' carriers 
comes from the Co and O Is core level positions. These 
measure the variation of the local electrostatic potential 
in the supercells. Supercell A has two different types of 
Co, the four Co coordinated by a Na in the prism im- 
mediately above or below and the two without. In the 
unrelaxed cell those with Na are 0.033 eV higher than 
the Co without. This variation is higher than the en- 
ergy difference between the e' band maximum and Ep. 
Secondly, while supercell B has all Co equivalently coor- 
dinated, the O sites differ in coordination in both super- 
cells. The variation in O core level position is 0.06 eV in 
both A and B. This variation may be significant because 
the bands are in fact hybridized Co-0 bands. The vari- 
ations in core level position are largely preserved in the 
relaxed cells, and structural inhomogeneities in the Co-0 
bond lengths and angles are introduced, which will also 
contribute to scattering. The variation of Co Is levels 
is reduced to 0.029 eV in supercell A, i.e. only slightly 
smaller than in the unrelaxed cell, and still larger than 
the e' band maximum relative to Ep, while the variation 
in the O Is positions increases to 0.07 eV.^* Thus the po- 
tential variations due to Na disorder are strong enough 
to localize the eg statesi2& 

The relaxation in supercell A is larger than in supercell 

B. This is because the Co and nearby Naco ions are both 
positive and move apart by 0.02 A, while the nearest 
neighbor Na-0 distance contracts by between 0.006 A to 
0.01 A, depending on the site. There are also distortions 
of the O cages. Measured by the difference between the 
shortest and longest Co-0 distance for a given Co site, 
these range from 0.01 A(2/3 of the Co in B) to 0.015 A 
(for the Co with Naco neighbors in A). 

In any case, the results show that scattering due to 
Na disorder is strong enough to localize the e' carriers 
at a; = 2/3. The actual position of the e' band maxi- 
mum relative to Ep is a function of doping, as is the Na 
distribution. The linear size of the e' Fermi surface, as 
predicted by band structure calculations varies from zero 
at high X to approximately 0.18% of the dimension of the 



zone at 2; ~ 0.3i^ At a;=0.3, the e' maximum lies 0.1 
eV above Ep- Observation of this Fermi surface would 
require a mean free path for e' carriers of at least the in- 
verse of this reciprocal space size. At x ~ 2/3 this would 
imply iegi > 50 A, where l^g' is the mean free path for 
the heavy carriers on the small Fermi surfaces. In met- 
als with Fermi surfaces having similar orbital character 
(here Co i2g), the constant scattering time approxima- 
tion is reasonable. Then i=vpT, where t is a scatter- 
ing time, vp is the Fermi velocity on a given section of 
Fermi surface and £ is the mean free path for carriers on 
that section. Taking the ratio of the velocities for the 
two sheets along the T-M line, the implication is that a 
mean free path lag ~ 250 A, for the main Fermi surface 
would be needed in order for the small section to be seen. 
Note that the small sections are elliptical in shape, elon- 
gated along T-M, so the criterion in the tangential direc- 
tion would be the same as the velocity is larger, and the 
size smaller, in proportion. Within kinetic transport the- 
ory, the conductivity is given by cr oc N[Ef)v'^t, where 
N{Ep) is the density of states at Ep and vp is the aver- 
age Fermi velocity in the direction in which the conduc- 
tivity is measured. Thus the large difference in the Fermi 
velocity between the large and small sections means that 
the conductivity is dominated by the large section, re- 
gardless of whether the small sections are present or not. 
Thus the measured resistivity is, for practical purposes, 
governed by the transport in the ag section. Depending 
on the doping level, Na2;Co02 samples of the type used 
in photoemission have been reported to have ratios of 
the resistivity at 300K to the residual resistivity in the 
range of 20 - SOfiiSi and resistivity saturation at or be- 
low ~ 600K. These values are inconsistent with a mean 
free path at low temperature as long as 250 A. We note 
that although the small sections become larger as x is 
lowered, they do so more slowly than in a rigid band pic- 
ture and the large ratio of the vp on the two surfaces is 



maintained. ^'^ 

Very recent high resolution photoemission workii 
shows dispersions of the Ug and e' bands below Ep that 
have a noticeable anti-crossing along T-K, for which the 
one electron e' and Ug bands have different symmetries 
and would not mix. The observed mixing is strong. Fur- 
thermore, a non-dispersive spectral weight is seen just 
below Ep in the region where the e' band is above Ep 
in band calculations. Both of these features are natural 
if scattering strong enough to localize the e' pockets is 
assumed. 

To summarize we find that disorder in the Na layer 
of Naa;Co02 produces sufficiently strong potential varia- 
tions to localize the e'g Fermi surface pockets at least for 
high X. Other scattering mechanisms and defects will en- 
hance this tendency. We argue that the non-observation 
of the small pockets and other features seen in recent pho- 
toemission experiments can be qualitatively understood 
within this framework, specifically that the carriers as- 
sociated with the small pocket are localized due to scat- 
tering related to Na and other disorder. In the hydrated 
superconducting compound, each Na is coordinated with 
four II2O molecules, is far from the C0O2 sheets and the 
Na is in a more ordered pattern with respect to the Co 
atoms. One may speculate that this leads to a sufficient 
reduction in the potential scattering due to the Na, and 
that this could lead to the re-appearance of the e' sec- 
tions. 
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